SOLAR  RADIATION  PRESSURE  ON 
HIGH  ALTITUDE  SATELLITE 


StOSHE  CARMEU 

(.7  SERAI  PHYSICS  RESEARCH  LABORATORY 


4 

9  -T  f# 

NATIONAL  TKHNICAl 
INFORMATION  SERVICE 

PROJECT  NO.  7114  *• 


D  D  c 


T 

Approved  for  public  release;  distribution  unlimited.  \_ 

ll 

)Tm  is 

c 

AIR  FORCE  SYSTEMS  COMMAND 

United  States  Air  Fares 


NOTICES 


When  Covcmmcnt  drawings,  specifications,  or  other  data  are  lived  for  any  purpose  other  than  in 
connection  with  a  definitely  related  Government  procurement  operation,  the  United  States  Government 
thereby  incurs  no  responsibility  nor  any  obligation  w  iiatsoevcr;  and  the  fact  that  the  Government  may 
have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is 
not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  iny  patented 
indention  that  may  in  any  way  be  related  thereto. 


Agencies  of  the  Department  of  Defense,  qualified  contractors  and  other 
government  agencies  may  obtain  copies  from  the 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 


This  document  has  been  released  to  the 

CLEARINGHOUSE 

U.  S.  Department  of  Commerce 

Springfield,  Virginia  22151 

for  sale  to  the  public. 


Copies  of  ARL  Technical  Documentary  Reports  should  not  be  returned  to  Aerospace  Research 
Laboratories  unless  return  is  required  by  security  considerations,  contractual  obligations  or  notices  on 
a  specified  document. 

AIR  FORCE!  8*5-72/200 


□  D  : 


1  *  £ 
;  I  1 


$  3 

ill 


£  o 


'35  * 
£  • 

B  5 
&  3 


cri 


ARL  72-0034 


SOLAR  RADIATION  PRESSURE  ON 
HIGH  ALTITUDE  SATELLITE 


MOSHE  C  ARM  ELI 

GENERAL  PHYSICS  RESEARCH  LABORATORY 


FEBRUARY  1972 


PROJECT  7114 


Approved  for  public  release;  distribution  unlimited. 


AEROSPACE  RESEARCH  l  ABORATOR1CS 
AIR  FORC  E  SYSTEMS  COMMAND 
UNITED  STAHS  AIR  FORCE 
VVRIGH'i  PVm  RSON  AIR  FORCE  BAS!  .  OHIO 


FOREWORD 


This  technical  report  presents  results  of  resear. a  carried  out 
by  Dr.  Moshe  Carmeli  of  the  Theoretical  Physics  Group  of  the  General 
Physics  Research  Laboratory,  Aerospace  Research  Labo  f.tories,  Project 
71L4.  The  research  was  initiated  by  the  suggestion  of  Dr.  Kenneth  E. 
Kissell,  Director  of  the  General  Physics  Research  . Moratory,  in  support 
of  the  Optical  Properties  of  Space  Objects  (OPOS)  program. 


II 


ABSTRACT 


In  this  report  an  exact  expression  to  the  solar  radiation  pressure 
on  a  prolate  spheroid  is  given. 


INTRODUCTION 


Recent  studies  by  Fea^,  and  by  Fea  and  Smith^,  have  shown  the 
existence  of  an  unexplained  acceleration  on  spacecrafts  which  are 
balloon  satellites  of  large  area  to  mass  ratio.  It  was  suggested  that 
the  unexplained  acceleration  might  be  caused  by  solar  radiation  pressure. 
There  is  evidence^ that  one  of  the  satellites  showing  this  anomalous 
acceleration  is  no  longer  spherical,  and  that  it  is  probably  shaped 
like  a  prolate  spheroid.  Previously,  it  had  been  assumed  that  the 
solar  radiation  scattered  by  the  satellite  is  symmetrical  about  the 
satellite-sun  line,  an  assumption  which  no  longer  holds  for  a  prolate 
spheroid.  It  is  therefore  expected  that  additional  perturbations  of 
the  orbit  will  arise. 

Consequently,  Smith  and  Fea^  developed  a  perturbation  method  to 
calculate  me  radiation  pressure  on  a  prolate  spheroid.  To  this  end 
two  major  assumptions  were  made:  (1)  the  effective  direction  of 
reflection  of  the  specular  flux  is  determined  by  Snell's  law  on  the 
incident  ray  that  passes  through  the  center  of  the  satellite;  and  (2) 
that  the  magnitude  of  the  flux  reflected  in  this  direction  approximates 
to  that  which  would  be  reflected  by  a  sphere  of  surface  area  equal  to 
that  of  the  spheroid.  Both  assumptions  seem  not  be  valid  in  general. 

In  this  paper  we  give  an  exact  expression  to  the  radiation  p  essure 
on  a  prolate  spheroid  without  making  any  one  of  the  mo  assumptions 


mentioned  above. 


PRELIMINARIES 


It  will  be  assumed  that  the  solar  radiation  in  the  vicinity  of 
the  satellite  is  homogeneous.  Also  it  will  be  assumed  (see  Introduction) 
that  the  satellite  has  the  shape  of  a  prolate  spheroid.  The  relative 
orientation  of  the  satellite  with  respect  to  the  direction  of  light 
is  then  determined  by  the  angle  between  the  semi-major  axis  and 
direction  of  light. 


We  will  use  Cartesian  coordinates  defined,  as  usual,  by  x  »  R 
sin  8  cos  cp,  y  ■  R  sin  8  sin  cp,  and  z  »  R  cos  8.  Let  the  prolate 
spheroid  be  located  at  the  origin  of  coordinates  and  described  b>  the 
equation 


x2  +  y2  z£ 

a2  c2 


1. 


(1) 


Also,  let  the  light  rays  be  parallel  to  the  x-z  plane,  having  the 
direction  given  by  (see  Figure  1) 


$  “  (-sin  0  ,  0,  -  cos  8  ) 

0  O' 


(2) 


radiation  pressure 

An  infinitesimal  light  beam  reflected  from  a  surface  element  dS 
of  the  spheroid  will  apply  a  force  given  by 

dF  =  21  t,n-N)  dS,  (3) 

-4 

where  l  is  the  light  density,  n  is  a  unit  vector  perpendicular  to  the 
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surface  element,  and  ft  is  the  direction  of  light,  which  is  given  by 
Eq.  (2).  The  geometry  of  the  surface  of  the  prolate  spheroid  determines 
both  n  and  dS.  They  are  given  by0 


n  *  a  g  ‘  sin  0  (c  sin  8  cos  cp,  c  sin  9  sin  cp,  a  "os  0), 


dft  »  g  ^  d6  dcpn, 


(4) 


(5) 


whtre 


2  2  2  2  2  2 
g  =  a '  sin  8  (a  cos  0  +  c  sin  0). 


(6) 


The  result  is 


.d,  21  a  sin  8  (c  sin  0  sin  0  cos  cp  +  a  cos  0  cos  9  ) 

or  “  "  _  o _ o 

,  2  2*,  2  .  2  _  v  ^ 

(a  cos  9  +  c  sin  0) 


x  (c  sin  G  cos  cp,  c  sin  8  sin  cp,  a  cos  0)  d0  dcp.  (7) 
The  total  force  F  is,  accordingly, 


F  - 


-♦ 

dF, 


(8) 


Here  a  stands  for  the  surface  snaded  by  light  on  which  the  integration 
should  be  carried  out. 


BOUNDARY  CONDITION 

The  boundary  condition  (integration  limits)  of  the  integral  of 
force  (8)  is  determined  by  the  curve  0  which  is  the  boundary  of  the 
o'  - A -  !  area  a.  This  curve  is  giv* n  by  the  equation 
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N*n(0,cp)  *  0. 


(9) 


Using  Eqs.  (2)  and  (4)  for  N  and  r.  one  obtains  explicitly  for  C  the 
equation 

c  sin  0  sin  0  cos  cp  +  a  cos  0  cos  0  *=  0.  (10) 

o  o 

This  is  the  equation  of  a  curve  given  in  terms  of  spherical  coordinates, 
cp  *°  cp  (0 ) ,  or 

cp  *  +  arc  cos  (-  ■§  cot  0q  cot  0).  (11) 

Notice  that  the  curve  cp  ■  cp ( 9 )  is  located  in  the  plane 

c  sin  0  x  +  a  cos  0  z  *  0.  (12) 

o  o 

In  Figure  2  we  give  the  boundary  limits  where  the  integration 
on  the  angles  0  and  cp  should  be  carried  out. 


TOTAL  FORCE 


The  force  integral  (8)  can  now  be  given  explicitly  as 


„  OT  rr  sin  0  (c  sin  sin  0  cos  cp  +  a  cos  0  cos  0)  cos  cp 

F  *  ’  21  ac  JJ  - - - 2__ - - - - - ^ - d6 


X 


f  2  2  ,  2  ,  2  . 

(a  cos  0  +  c  sin  0) 


(13a) 
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„  ffsin  0  (c  sin  0  sin  0  cos  cp  +  a  cos  0  cos  0)  sin  cp 

F  »  -21  ac  JJ _ _ _ o _ _ _ o _ _ _ 


(a  cos  0  +  c  sin"  0)^ 


dcp  d0 

(13b) 


Fz  *  "  21  *  JJ 


sin  0  cos  0  (c  sin  0  sin  0  cos  cp  +  a  cos  9  cos  0) 

o  o 


,2  2  2  .  2  .  ^ 

cos  8  +  c  sm  0) 


dq-  dB 


(13c) 
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As  can  be  seen  from  Figure  2.  the  limits  of  the  double  integrals  in 

Eqs.  (13)  should  be  taken  from  cp  =  -  rt  to  cp  =  +  rr  for  Os  6  s  8^,  and 

from  cp  »  -  arc  cos  (-—  cot  0  cot  0)  to  »  *  +  arc  cos  (-5  cot  0  cot  0) 

co  co 

for  Q^sesrr-©^.  The  angle  0^  is  that  0  for  which  the  angle  cp  of 
the  curve  C  satisfies  w(0^)  -  rr. 

The  force  term  F  ,  Eq.  (13b),  can  easily  be  shown  to  be  equal  to 

zero. 


CASE  OF  A  SPHERE 


The  force  terms  given  by  (13)  have  a  particularly  simple  form 
in  the  case  of  a  sphere  (i.e.  a  *  c).  One  obtains  in  this  case 


F  **  ■  -  21  a^  ffsin  0  (siri  9  sin  0  cos  cp  +  cos  9  cos  0)  cos  cp  dcp  d0 
x  J  J  o  o 

(14a) 


F  S  -  0 

y 


(14b) 


S  2  r  (* 

*  -  21  a  JJsin  9  cos  0  (sin  sin  9  cos  cp  +  cos  9o  cos  0)  dcp  dO 

(14c) 


One  can  simplify  these  integrals  by  going  into  new  variables  0 ' , 
c;'  m  which  the  integration  limits  are  not  dependents  on  each  other. 
To  this  end  one  rotates  the  coordinate  system  around  the  y-axio: 

x  «  x'  cos  9  +  z  '  *»  Ln  9 

■j  o 

z  "  -x '  sin  0  v  2'  cos  -j  .  (11)) 

o 


A  simple  alcciatlon.  using 


the  relations  x  =  a.  sin  9  cos  cp,  y 


a  sin  0 


sin  cp,  z  *  a  cos  9  and  x'  =  a  sin  0'  cos  cp',  y'  =  a  sin  9'  sin  cp ' , 
z'  «  a  cos  9',  then  shows  that  0  and  cp  are  related  to  the  new  variables 
0 '  and  cp '  by 

sin  8 'cos  cp'cos  6  +  cos  8 'sir  9  =  sin  0  cos  cp, 

o  o 

-sin  8'  cos  cp'  sin  0  +  cos  8'  cos  6  «  cos  0,  (16) 

T  o  o  ’ 

and 

sin  8 'd0  '  dcp'  a  sin  8  d0  dcp.  (17) 


As  a  result  one  obtains  for  the  force  components  acting  on  the 
sphere 

TT 

S  V  r"?  f4-n 

*X  V  ^  &  J  j  (sin  C0S  Cp#  COS  60  +  C0S  0/  sl-n  0Q)  COS  0/ 

9=0  Cp=-TT 

sin  0  '  dcp '  d0  ' 


-21  sl~  ,7  r  (-sin  9'cos  cp  sin  8  +  cos  8'  cos  0  )  cos  8' 

J  J  o  o 

9  =  0  Cp“-TT  ,  . 

*  sine'  dcp'  d0 

(18) 


The  result  is 


F  ^  =  -(4/3)tt  a'  I  sin  0 
x  o 

F  S  =  0 

y 

F„^  =  -(4/3 )n  a~  I  cos  e 
L  o 


(19) 
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Figure  1.  Prolate  ellipsoid  and  light  rays  in  the  x-z  plane  along 
the  vector  ft  =  (-sin  8  ,  0,  -  cos  8  ).  The  curve  C  is 

•4-4  ®  0-4 

defined  by  N*n(8  ,cp)«*0,  where  the  vector  n  is  normal  to 
the  surface  of  the  ellipsoid. 


Figure 


4> 


2.  Boundary  limit  diagram.  The  boundary  of  integration  is 
given  by  cp  *>  +  n  for  U  *  9  &  0  and  cp  ( 9 )  «  Hr  arc 
cos(-|.  cot  6q  cot  8)  for  8^  *:  0  s  tt-0^  The  angle  01  is 


defined  by  8,  »  arc  tg  cot  0  ). 
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